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Abstract. 

We have studied the conformal models WD^\ n = 3,4,5,..., in the presence of 
disorder which couples to the energy operator of the model. In the limit of p ^> 1, where 
p is the corresponding minimal model index, the problem could be analyzed by means of 
the perturbative renormalization group, with e - expansion in e = -. We have found that 
the disorder makes to flow the model WD%) to the model WD^ -1 ^ without disorder. 
In the related problem of N coupled regular WD^ models (no disorder), coupled by 
their energy operators, we find a flow to the fixed point of N decoupled WD&~ 1 '. But 
in addition we find in this case two new fixed points which could be reached by a fine 
tuning of the initial values of the couplings. The corresponding critical theories realize 
the permutational symmetry in a non-trivial way, like this is known to be the case 
for coupled Potts models, and they could not be identified with the presently known 
conformal models. 
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1 Introduction 



With the present knowledge of the effect of disorder on the critical behavior of statistical 
systems it could be fair to say that any extra theoretical model will be interesting in 
which these effects are nontrivial and can be studied reliably 

In this respect we find the set of conformal theory models WD^\ n = 3,4,5,..., 
to be interesting. Here p is the index of the corresponding unitary series of minimal 
models. These models are based on the corresponding W - chiral algebras. They have 
been defined in the paper Q. 

The simplest model in this series WD^ will be sufficient to give a preliminary dis- 
cussion of our approach and results. 

The model itself, which has not much been used in statistical physics applications, 
will be described in the next Section. For the moment we shall observe that a particular 
primary operator in this theory, which is naturally identified with the energy operator, 
e(x), possesses conformal dimension 

A, = i (1.1) 

in the limiting theory WD^\p — > oo. Physical dimension of e will be twice bigger, 
A p e h = 2A £ = 1. For p finite, A e is of the form: 



where 



If we couple disorder to e(x) 



p+l 



[1.3) 



d 2 x m(x)e(x) (1-4) 



with m(x) being random valued, then the disorder produced interaction of the replicated 
theory: 

N 

g / d 2 xJ2 £ a(x)e b (x) (1.5) 
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will be marginal in the limit p — > oo, and will be slightly relevant in the model WD^ 
with p large but finite. In this case the effect of disorder can reliably studied by the 
perturbative renormalization group, by developing in e. 

We observe also that the problem formulated in this way will be defined better, 
theoretically, compared to the much used Potts model, or the minimal conformal theory 
models, with disorder [2,3]. In this latter case the disorder is marginal in the case of the 
Ising model [4] , and the e parameter of the Potts model, in which, similarly, the disorder 
is slightly relevant, corresponds to the deviation of the central charge, or of the number 
of components, from the Ising case. The point of difference is that in this case there 
are no unitary models close to Ising which would correspond to infinitesimal values of 
e. Ising model is not the accumulation point of the unitary series. Next to Ising, the 
3-states Potts model, lies at numerically small but finite distance in e. For this reason, 
the perturbative expansion in e is likely to give numerical results, in case of Potts model 
with disorder, not analytic ones. 

This is different in case of the model WD^ where the disorder is marginal at the 
accumulation point of unitary series, p = oo. In this case for any small value of e ~ 1/p 
there would be a physical, unitary model behind. In this sense, e will be a well defined 
analytic parameter. One could expect that the perturbative expansion in e would provide 
in this case an analytic expansion in e = l/(p+l) of unknown yet exact functions (central 
charge, dimensions of operators). Like this is actually the case for the minimal model 
perturbed with the operator 3i i [5,6]. For these reasons we consider the model to be 
interesting. 

Another important point of difference with the Potts model is that in the operator 
product decomposition of e(x)e(x') there appears an extra slightly relevant operator, we 
shall note it <f>(x). The operator algebra has the form: 

e(x)e(x>) = j—^il + ...) 

+ | X _^W ^ + -) + \ x - jji-A, + -) ( L6 ) 

The detailed information on this decomposition shall be given in the sections 3. For 



the moment we shall observe only that one of the operators in (1.6), let's say (f>'(x), 
is irrelevant, so that it can be dropped in the RG calculation. While another, <j)(x), 
is slightly relevant = 1 — Ae. This requires, in the RG calculation, to add this 
extra operator to the perturbative action, so that the replicated theory will have the 
perturbation of the form: 



with the initial values, in the RG sense, g = —go, go > 0, which is due to disorder, and 
A = Ao = 0, initially zero, but produced further spontaneously by the interactions. 

As will be shown in the next Sections, the theory with A pert in (1.7) has a stable fixed 
point at 



and it corresponds to N decoupled models WD$ , so that, in the limit N = which 
corresponds to problem with disorder, one gets a flow: 



One gets this flow for N = 0, in the coefficients of the (3 - function, and for the initial 
conditions g < 0, A = 0. There are no other fixed points in this case (N = 0) in 
addition to (1.8) and the trivial fixed point, unstable, g — A = 0. Yet in the related 
problem of N coupled regular models (no disorder), for N = 3 or bigger, the (3- function 
contains, in addition to (1.8), two other fixed points, both stable in one direction and 
instable in the other. They correspond to multicriticality. To reach them one has to 
have both initial couplings g , A ^ 0, and, in addition their ratio have to be fine 
tuned. These extra fixed points are nontrivial, in the sense of a nontrivial realization at 
these points of a permutational symmetry. They are also nontrivial in the sense that the 
corresponding couplings (g*~\ A^), (gi 2 \\^) are all non-vanishing. We expect that 
these points would not be identified with any presently known conformal field theory, like 
this appears to be the case for the coupled Potts models [2,7]. As has been argued above, 
we expect that in this case, of coupled WD^ and more generally WD^ models, our 
calculations provide first analytic corrections in e ~ ^ to the development of yet unknown 
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(1.7) 



0, g* = 



(1.8) 



WD^with disorder — > W^ p ^regular, without disorder 



(1.9) 



exact functions of e, or p. In case of the Potts models, the corresponding corrections in 
e have to be considered as only numerical. 

The rest of the paper is organized as follows. In Section 2 we introduce the Coulomb- 
Gas representation of the WD^ model. Then, in Section 3, we study the effects of weak 
bond randomness on this model and the related problem of N coupled WD^ models. 
The renormalization of the couplings and of the energy operators are computed at one 
loop order. In Section 4 we study the behavior of the RG-flow in the case of disordered 
system and of N coupled systems; we show in particular the presence of non trivial fixed 
points in both cases. We give the analytic expansion in e of the central charge and of the 
dimension of energy operators at these points. In Section 5 all the results we have found 
for the WD { i ] model are generalized for the whole family of WD&> models. We discuss 
what has been obtained in the Conclusions. 



2 Coulomb gas of the model WD% 

To make our presentation self-contained, and also because the models {WD^} have not 
much been used so far in the statistical physics applications, we shall reproduce in this 
Section the results of the paper [1], by giving the details on a particular model of our 
interest, WD^K 

This model could be realized by a 3-component Coulomb gas, with the stress-energy 
operator T(z) taking the form: 

T(z) = ~ : d(p(z)d(p(z) : +ia o d 2 (2.1) 

Here <p(z) = {<pi(z), (p2(z), <P3(z)} is a set of three free fields put into a vector, with the 
correlation functions normalized as: 

< Lp a (z,z)Lp b (z',Z>) >= 2 log | l —^5ab 

\z — z y 

= (2 log — ^ + 2 log ^-^)5 ab (2.2) 
z - z' z — z' 

In the following, and in the above formula (2.1), we often suppress the dependance of 
ip a (z,z), on z, which will be implicit. 

5 



The vector a in (2.1) corresponds to a presence of the background charge operator, 
V-2a {R)-> putted at infinity, R — > oo [8]. 

In addition to T(z), the model contains two other operators in its chiral algebra: 
W 3 (z), W 4 (z), with conformal dimensions 3 and 4. They could also be expressed in 
terms of polynomials in derivatives of fields {(p a (z)} [1]. The expressions are longer, 
but they will not actually be needed in the present analysis. Important is that W 3 (z) 
and W^z) exist and classify, together with T(z), all the fields (operators) of the model 
in terms of primaries and descendants of the chiral algebra. Then the usual methods, 
combined with the available Coulomb gas representation, define dimensions of primary 
operators (Kac formula) and their correlation functions. In this sense the conformal 
theories of WD$ are fully defined [1]. 

In the Coulomb gas of the model WD^ there are 3 screening operators "+" , and 3 
screening operators "-" : 

V-{z) = V sz {z) = e^^ (2.4) 
"a = a±4, (4) 2 = 1 (2.5) 

The "length" of the screening operator vectors, ol±, are fixed by the condition that the 
conformal dimension of V^(z) is equal to 1. In general, for a vertex operator 

V s (z) = e iS ? (2.6) 

its conformal dimension with respect to the stress-energy operator T(z), eq.(2.1), is given 
by: 

A d = A(Va) = a 2 - 2aa = a 2 - 2aa cos 9 (2.7) 

Here 6 is the angle between a and a. For the geometry of screenings of WD^ the 
angle between {a a } and a is the same, Fig.l. One gets, as a condition A(V r a ± ) = 1, 

a 2 a — 2a a a cos = 1 (2.8) 

( a a)l,2 = Oi± — ttQCOS© ± \J Oq COS + 1 (2.9) 
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The unit vectors {e a } shown in the Fig.l correspond to simple roots of the classical 
Lie algebra D 3 . As it was said above, the vector a is supposed to be "equally distant" 
from e*i, e*2, e*3. Assuming this condition, for which we will not go into details here, cos0 
will be defined by the geometry in Fig.l. One finds: 

cos0 = ^= (2.10) 



By the formula (2.9), the lengths ct± of the screenings in eq.(2.5) will be functions of 
a only. This implies that in the model WD^ the orientational geometry is fixed, but 
there is one free parameter in the lengths of the vectors: the screenings {af} and the 
background charge c?o- 

The primary operators of the model are represented by the vertex operators 

V {z) = e l ^ z) (2.11) 

The allowed values of the vectors (5 are defined by the degeneracy condition of the modules 
of Vg{z) with respect to the chiral algebra. They are found to be given by [1]: 

-> -> 3 1 — n' 1 — n 

(3 = /V^iX^WK.ns) = XX — T JLa ~ H o~ ~ a +^ ( 2 - 12 ) 

a=l ~ ~ 

Here {uj a } is a set of three vectors which are dual to the unit vectors of the screenings 

(uj a ,e b ) = 5 ab (2.13) 

Using the general formula for the dimensions of vertex operators, eq.(2.7), one then gets 
from (2.12) the Kac formula of the model which defines the set of dimensions of the 
primary operators {<P{ n ' 1 ,n 1 ){n' 2 ,n 2 )(n' 3 ,n z ){.z)} : 

K ,nO(«> 2 )(n>3)W « % i , ni)K ,„ 2)(n ,, n3) (^) (2-14) 

A(„' iini )(„^ n2 )( n ^ n3 ) = - 2a /3(...) 
= E(««. n a )f^f + 2 £ ««, n>K, n b )^p± 

a 4 a<b 4 

+ n a )^^ (2.15) 



Here 

u ( n ' a i n a) = (1 - n' a )a- + (1 - n a )a + (2.16) 
In (15) we have used in addition the decomposition: 

2a = K + a-)£<4 (2.17) 

a 

which is easily verified by multiplying both sides with a vector e h and using: 

2cJ efe = 2a cos = a + + a_ (2-18) 



- by eq.(2.9), and (uj a eb) = S a b, eq.(2.13). The unit vectors {e a }, which are defined by 
the Fig.l, when expressed in components will be given by: 



Si = ^(0,-1,1), e 2 = -L(1,1,0), e 3 = -L(-l,l,2) 
The dual vectors {uj a }, defined by the eq.(2.13), are found to be equal to: 

^ = ^(0,0,1), £3 2 = -^=(1,1,1), U 3 = ^=(-1,1,1) 



Then one finds: 



(w ,t3 6 ) 



' 2 1 ^ 

1 2 I 

2 2 

i 1 3 

\ 1 2 2 / 

)(P) 



Finally, the Kac formula of the model WD 3 takes the following form: 
A 



1' ' ' ^2 i 3/ / , /n.n.2 , 3, , , ,>> 2 
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+ -w(n / 1 ,n 1 )w(n 2 ,n 2 ) + -w(ni, ni)w(n 3 , n 3 ) + -u(n' 2 , n 2 )u(n' 3 , n 3 ) 

3 3 
-(«+ + a^)(2u(n' 1 ,n 1 ) + -u(n 2 ,n 2 ) + -u(n 3 ,n 3 )) 



u(n' a) n a ) are defined by eq.(2.16). One checks that 

^2(u a ,uj b ) = (4,3,3) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



which we have used in the last term in (2.22). 



As it was mentioned above, the model contains one free parameter. This will be either 
a , or a + , or ck_, or the central charge of the model, which, for T(z) in eq.(2.1), will be 
given by: 

c = 3 - 24a 2 (2.24) 

One observes that the lengths of the screenings "+" and the screenings "-" , a + and a_ 
in eq.(2.9), are related: 

a + + a_ = 2a cos @ (2.25) 

a+a ^ = -1 (2.26) 

Similar to the case of the basic conformal theory of [9], one gets a unitary set of models 
if the parameter (a+) 2 is restricted to the discrete values: 



K) 2 = ^ (2.27) 



By the relation (2.26) one will have also 



(«-) 2 = ^7 ( 2 - 2 8) 

For WD^\ p = 5,6,... [1]. In the following we shall restrict ourselves to this unitary 
series. This explains for the extra index p of WD^\ 

We shall next be interested in this model in the limit of p large, p> 1. First, when 
we take p — * oo, the Kac formula (2.22) simplifies and takes the form: 

a+ = l, a = -l (2.29) 

u{n' a ,n a ) =n' a -n a (2.30) 

13 3 

A (n>1 ) K , n2 )( n>3 ) = -K - m) 2 + -K - n if + g K - n s) 2 

- n i)( n 2 - n 2) + - - n s) + - n 2 )(ng - n 3 ) (2.31) 

In general, not specifically for p large or infinite, the model contains a symmetry, in the 
set of its primary operators and their dimensions, w.r.t. permutation of indexes 2 and 
3. This corresponds to the reflectional Z 2 symmetry of the geometry of the screening 
vectors {e a } in Fig.l. 
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In the set of dimensions A( n / )ni )( n ^ n2 )( n / jri3 ) there is a subset which is not Z 2 degen- 
erate and which corresponds to the singlet operators. For instance, singlet operators are 
{0(n' 1 ,ni)(i,i)(i,i)}- The lowest dimension operator in the singlet sector, next to the identity 
operator I = 0(i,i)(i,i)(i,i)(A (1)1)(1)1)(1)1) = 0, by eq.(2.22)), will be the operator 

0(2,i)(i,i)(i,i)Or) (2.32) 

It is natural to identify it with the energy operator of the model, 

e(x) = <i>(2m,i)(i,i)(x) (2.33) 

If we look now at the limiting case of p — > oo, we shall find, by the formula (2.31), that 

A £ (p^oo) = i (2.34) 

For p finite, by the formula (2.22), one finds: 

A £ = A (2>1)(1>1)(1>1) = - 2 (2.35) 

Substituting, by eq.(2.28), 

a 2 = = 1 — = 1 - e (2.36) 

p+1 p + 1 



where we have defined 



one obtains: 



1 



p+l 



(2.37) 



A £ = \ - \e (2.38) 

In our present study we are going to couple disorder to the energy operator e(x), defined 
in eq.(2.33). In the related problem of coupled regular models, we shall couple the models 
between themselves by their energy operators: 

N 

Apert = g^2 d 2 x e a (x)e b (x) (2.39) 

As has been discussed in the introduction, in both cases one gets a problem with a slightly 
relevant perturbation, if A e is given by the eq.(2.38) with e small, or p large, eq.(2.37). 
These problems will next be studied by the methods of the perturbative RG. 
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3 1-loop RG equations for WD$ 



3.1 Renormalization of the couplings 

We initially consider N regular WD^ models coupled by the energy-energy interaction 
(2.39): 

N N . 

A = Y / 4 a) +goY, d 2 xe a (x)e b (x) (3.1) 

Aq 1 ^ being the conformal action corresponding to a sing le WDf model. The action (3.1) 
describes a conformal field theory perturbed by a slightly relevant term quadratic in the 
energy operator; such a problem can be reliably studied by means of perturbative RG 
with e-expansion. The RG scheme requires that all the relevant terms produced by the 
energy-energy interaction have to be added to (3.1) and the algebra of the enlarged set 
of perturbing fields have not to present other relevant operators. As shown below, the 
O.P.E. of e(x)e(x') contains, apart from identity, one and only one (slightly) relevant 
operator, namely 0(i,i)(2,i)(2,i) = <f>, whose dimension is A( 11 )( 2i i)(2,i) = 1 — 4e. This 
implies that in the perturbative computation there will be diagrams (generated by terms 
in (2.39) with the same replica index) which produce, apart from trivial or irrelevant 
contributions, the term: 

N 

d 2 x<t> a {x) (3.2) 

a J 

We therefore consider the more general action: 

N N N 

A = Y / A i a) +goY, d 2 x e a {x)e h {x) + A £ / d 2 x <f> a {x) (3.3) 

The problem of a single WD^ model perturbed by has been considered in [1], where 
the existence of a non trivial infrared fixed point to which the system flows has been 
shown. The conformal field theory associated to this point corresponds to the WD^~^ 
model. 

In order to investigate the energy algebra, we study the four point correlation func- 
tion G(x) = {e(Q)e(x)e(l)e(oo)) , which can be decomposed in a sum over s-channel 
diagrams corresponding to insertions of different operators. Introducing the Coulomb 
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— # — # 

gas representation (c/r. section 2) for G(x) and defining for simplicity /3(2,i)(i,i)(i,i) = Pe, 

2a - /3(2,1)(1,1)(1,1) = Af, /9(1,1)(2,1)(2,1) = aIld 2 «0 - /%,1)(2,1)(2,1) = we can write : 

G{x) oc f.J{V£(p)V A (x)Vz(l)V&(oo)Vr (/.OyfM^-^)^"^))^ 2 ^ 2 ^ 

(3.4) 

where the screenings V^ - , defined in (2.3), are integrated over the 2D plane and their 
number is determined by the charge neutrality condition : 

2(3 e + a_ (miei + m 2 e2 + m 3 e 3 ) = (3.5) 

which is satisfied for mi = 2, m 2 = 1, m 3 = 1. 

Each leading term in x — > limit originates in the integral representation (3.4) when 
a particular number g a of screenings (by eq.(3.5) q\ < 2, q 2 < 1, #3 < 1) approaches 
the origin. Since in this limit 

a 

the intermediate state vertex V 2 ^ +Q ,_ ^ QaSa corresponds to a primary field which belongs 
to the energy algebra. In particular we have noticed the presence of the vertex V 2 e +a-e^ 
which corresponds to the primary operator 0, according to 

2& + a_<r 1 = -a_^ + ^ (3.7) 

In appendix A we shall show in detail that is the only relevant operator in the energy 
product decomposition and that the algebra of and e, apart from irrelevant terms, is: 

I 

£{XH0) = 1^ + N 4A £ -2A^ + • • • 

<f>(x)e(0) = -^-0 + -..; (3.8) 

Df e and are the associated structure constants of the WD^ model. 

At one loop order the main RG characteristics of action (3.3), e.g. /3-functions or 
renormalization of operators, are easily obtained from eq.(3.8) (see appendix C). The 
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renormalized coupling constants g{r) and A(r) (r is the short distance cut-off) have been 
computed, and the correspondent /9-functions take the form: 

(2 - AA £ )g — (N — 2)g 2 - AB&g 

N — 1 

(2 - 2A,)A - -f\" - -^-Dtg 2 (3.9) 

where the coupling constants have been redefined as g — > g/An and A — > A/27T. In order 
to study the coupling flow induced by eq.(3.9) (section 4) the structure constants have 
been determined (see appendix B): 

^ £ = yf + 0(e), D*+ = -*= + 0( e ) (3.10) 

In section 4 we will investigate the properties of equations (3.9) for N = (disorder) 
and N > 2 (coupled systems). 



9 dln(r) 

, = 9A(r) 
dln{r) 



3.2 Renormalization of the energy operator 



In order to study the effect of the perturbation on the energy operators, we need to 
compute the renormalized operators e' a , which are expressed via the iVxiV matrix [Z £ ] ab 
by: 



£ 'a - ^[^e]o6 £b 



(3.11) 



As usual, we proceed perturbatively: computing contributions from each coupling term 
and rewriting bare quantities in term of renormalized ones, we have determined up to the 
first order the matrix 7^ = dln[Z] a b / dln(r) (see appendix C). The new critical exponents 
at the fixed points g*, A* of the RG flow will be given by the eigenvalues of the dimension 
matrix A ab = A £ 5 a b — labi.9*, A*) which takes the form: 



A 



ab 



~ 1 



9* 
9* 



9* 

A £2 + D*K 



9* 
9* 

+ D*K 



9* 



9* 
9* 
9* 



A £JV + D*K 



(3.12) 
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It is straightforward to see that the symmetric combination e s = S\ + e 2 + . . . + £jv and 
the antisymmetric ones e ai = S\ — e 2 , e a2 = e 2 — £3 ■ ■ ■ are eigenvectors of (3.12) and 
their dimensions turn out to be: 

A £s (g*,K) = A £ + £> £ %A* + (N - l)g, 

A eoi (^,A,) = A £a2 =... = Ac + D^-g* (3.13) 

In the next section we shall show that equations (3.9) admit fixed points with g* ^ at 
which N identical systems remain coupled; we expect therefore that the correspondent 
conformal field theory realizes in a non trivial way the permutation symmetry of iV 
identical objects. The splitting of dimensions in the energy sector A £s ^ A £a indicates 
that this is indeed the case. Although a similar result was obtained for iV coupled Potts 
models, in this case we have access to the analytic expansion in e of the dimensions of 
the energy operators and of the central charge (see next section) of new conformal field 
theories. 



4 RG-flow and central charge 

4.1 Fixed point structure and couplings flow 

The first step in the study of RG-flow is to find all points g*, A* such that (3 g (g*, K) = 

Px(g*,K) = 0. 

The equations (3.9) have been obtained for a generic number iV of coupled models. We 
consider firstly the quenched case which is obtained in the iV — > limit. The /3-functions, 
defined in (3.9), take in this limit the form: 

ft = 10e 9 + 2g 2 - 2i/|a 9 

The RG flow has, apart from the trivial solution g — A = 0, one stationary point: 

A* = 2Vl5e, g m = (4.2) 
14 



which is stable in all directions. This point has already been found in [1]: the associated 
critical theory is described by a pure WD^~^ model. Taking as initial conditions Ao = 
and go < 0, and assuming that higher loop corrections will not change the qualitative 
behavior of the flow near the two fixed points, the system flows toward point ( |4.2|) as 
supported by the numerical calculations shown in Fig. 2. 

Initial conditions A = and g < correspond exactly to the disorder produced cou- 
plings; so, the addition of a small bond randomness will drive the model WD^ to the 
model WD^ without disorder. 

We consider now N > 2, i.e. the case of coupled models: by (3.9) and (3.10) the 
/3-functions are: 

ft, = 10£ 9 - (N - 2)g 2 - 2^\g 

* = * A -^ A '-^i/!k (4 - 3) 

In addition to point (4.2), which remains a point of attraction of the RG-flow since the 
respective RG-eigenvalues don't depend on N, we have two new fixed points: 

The stationary points ( |4.4| ) and ([4.5|) constitute node points of the flow and they can 

be reached only by a fine tuning of the initial conditions. Studying in detail the flow 

diagram of (4.3) (shown for N = 3 in Fig. 3) we see that with the initial conditions Ao = 

and go ^ the system will flow far from our fixed points toward either a massive theory 

or another fixed point which cannot be seen at this order in perturbation theoryQ On the 

*It could be mentioned that a particular case of N — 2 coupled WD^ theories, with go ^ 0, Ao = 0, 
have been considered previously in the paper [11], with a conclusion that the theory is integrable. We 
only want to mention again that according to our RG analysis the second coupling, A ^ 0, will have to 
be admitted eventually into the action (see the corresponding flows in Fig. 3). 

This is in contrast with two coupled Virasoro algebra minimal models M p which were shown to be 
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other hand, in the region A > and of order e, the situation is rather different. There are 
two solutions of eq.(4.3), say <?+(A) > and g~(X) < 0, which are attracted respectively 
by the node point (4.4) and (4.5); if we take Ao > and g-(Xo) < 9o < g+{Xo), the system 
will flow toward the stable point (4.2). In this case the infrared limit of the system will 
be described by N decoupled WD^~^ models. 

It's important to note that, although quite similar, the flow is not symmetric along 
the A axis, and this will be at the origin of the difference between the values of the central 
charge at the two fixed points. This asymmetry can be explained by a simple physical 
argument. Indeed, when g > 0, the coupling of different models is always frustrated, 
while for g < the iV-system can arrange itself in a kind of ferromagnetic configuration 
in order to minimize the energy. For N = 2 (i.e no frustration) we recover the symmetry 
g — > —g, and the central charge has the same value at the two new fixed points (see next 
section). 

Inserting the values of the structure constants (3.10) in eq. (3.13), we have at the 
fixed points (^* 1 \a1 1 ^) and (gi 2 \\^): 

A (D,(2) = A £ + 5 (i ±N l 6 ) e 
\ V N - 1 + 6N 2 ) 

A (i),(2) = A e + 5|1 T W 6 7 \e (4.6) 
^ V N - 1 + 6N 2 J K ' 

where A« = A Es a (gi 1 \ X^) and Ag. = A^f, Al 2) ). By eq. (4.6) we have A« = 
Ai 2 ) and A« = A?> for all N. This result can be easilv explained for N = 2: in this 
case the two models are equivalent under the replacements g — > —g and E\ — > —E\. On 
the other hand, for iV > 2 no reason can be given, as it can be seen from the multiplicity 
of the antisymmetric energy combinations. We believe therefore that these equalities are 
accidental and they originate from the simplicity of the one-loop order computation, 
integrable [12] and studied in detail in [11]. In this last case the original action, with a coupling term g 
only, is stable with respect to RG evolution. 

Then, if the coupling A ^ has to be taken into the action, for the two coupled WD^ theories, 
eq.(3.3), the natural question would be if the presence of the second perturbative term will modify 
essentially the analysis of integrability. 
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4.2 Central charge 



In the previous section we have shown that the RG-flow of the quenched system (N = 0) 
and of the coupled systems (N > 2) exhibits non-trivial fixed points. A simple way 
for computing the central charge of the associated conformal theories is given by the 
Zamolodchikov's c-theorem || ; the theorem provides us with a function of the couplings 
c(g, A), to be defined below, whose value c(g*, A*) at the fixed point A*) is equal to 
the central charge of the corresponding critical theory. 

We define as Q(x) the trace of the stress-energy tensor. It's well known that 6(x), 
which is zero at the fixed point, is proportional to the perturbing terms of the theory. In 
our case we have from the action (3.3): 



4 



= f E £a(x)e b (x) + ^ E <j> c {x) (4.7) 



where the renormalization of the couplings g — * g/47r and A — > \/2ir is taken into account. 
The function c(g, A) is then completely determined by the following equations ||: 

P 9 9 -^f + fc^ = -24 < 6(0)0(1) > 

c(0, 0) = c pure (4.8) 



with Cp Ure and < 6(0)6(1) > respectively the central charge and the 6 field two-point cor- 
relation function of the unperturbed theory. Inserting eq.(4.7) and eq.(3.9) into eq.(4.8) 
and using the following relations: 

N N 

EE< e«(l)e6(l)e c (0)e d (0) >= 2N(N - 1) 

N N 

EE<^(o)^(o)0 Q (i)>=o 

c^d a=l 
N N 

EE<^( 1 )^(°)>=^ ( 4 - 9 ) 

a=l 6=1 

it's straightforward to verify that the function 

c(g, A) = c pure -N ^ W + 3A,eA 2 - {N - X M - 2) g 3 - ^f\ 3 - 5(JV - 1)Z^ 2 A 

(4.10) 
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satisfies eq.(4.8). The function (4.10) has been obtained when a number iV of WD$ ' 
models is considered and so c pure = Nc(WD^), where c(WD^) is the central charge of 
a single WD^ model, calculated in [1]: 

^Hsfi-^) (4 . n) 

In order to compute the central charge in the case of the disordered problem, we normalize 
the function c(g, A) by N and then we take the limit N — > 0; in fact, in terms of replicas 
this is exactly the limit in which the related quenched free energy is obtained. Using 
the values of the structure constants (3.10), the central charge Cdi S . at the infrared stable 
point (4.2) turns out to be: 

c dis . = Mm C( °' 2 f^ t] = c(WD<?>) - 1206 3 » ^WD^) (4.12) 

This result is consistent with what has already been said in the previous section: the 
infra-red behavior of a WD^ model with disorder coupled to the energy operator is 
described by a WD^~^ model. 

In the related problem of iV coupled models, the central charge c coup i, at the fixed 
point (4.2) is c coupL = c(0,2y/l5e) = Nc(WD^~^), i.e. the correspondent critical theory 
is described by N decoupled WD^~^ models. 

Finally we have access to the analytic expansion up to the third order in e of the 
central charges C\ and c 2 at the new fixed points (gi 1 ^, A^) and (gi 2 \ X^): 

ci, 2 ^ NciWD^) - N 

The eq.(4.13) represents the first analytic result for a new series of conformal theories 
which in addition to the W^-symmetry present a non-trivial representation of the permu- 
tation symmetry. 



60 I 1 =f (iV — 2) 1 



/ 6 
N - 1 + 6iV 2 



(4.13) 



5 Analysis of the WD)f> models 

We shall show that the results we have obtained in the previous sections still hold in 
the case of WD^ models. The construction of the Coulomb-Gas representation of these 
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theories is a direct generalization of what has been presented in Section 2. The WD^> 
model presents, in addition to the conformal symmetry, a series of additional symmetries 
which are generated by a series of local currents {W2k(%)}k=i,...,n-i and W n (x) with 
dimension = 2fc and A n = n respectively. It can be represented by an n-component 
Coulomb Gas. The (2.3), (2.4) and (2.5) still define the screenings "+" and with 
unit vectors e a which lie on a n-dimensional space and correspond to simple roots of Lie 
algebra D n . The primary operators are represented by vertex operators with 

P = P(n' 1 ,n 1 )(n' 2 ,n 2 )..(n' n ,n n ) = }^ ( g "~ H 9 a +' Wa 

a=l " ~ 

where the set of dual vectors {u) a }, defined by (2.13), have the quadratic form matrix 

F a b = (<^a> ^4) : 

F ab = 2a, a < b < n - 1; F an ^ x = F an = a, a < n - 1; 

F F - —■ F — U ~ 2 ^ 

r nn r n— In— 1 ^ r n—lm ^ v." / 

The central charge c(WD^) of the WD^ model has been calculated in [1]: 

Using (2.15) and ( |5.2D the dimension of primary operators can be easily calculated. In 
particular the dimension of the operator 3>(2,i),(i,i),..,(i,i)> naturally identified with the 
energy operator e (Section 2), is: 

1 2n- 1 ,„ , N 

A £ = A(2 ) i),(i,i),...,(i,i) = 2 — e ( 5 - 4 ) 

Therefore the disorder induced interaction (~ ee) is slightly relevant: we can study the 
WDffi model with disorder coupled to the energy operator and the related problem of 
iV-coupled systems using the same technique as the one we have exploited in the case of 
a WD^ ' model. The analogy with this case goes further: we show in appendix 1 that 
the operator </>(i,i)(2,i)(i,i)..(i,i) = with dimension 

= A(i,i),( 2 ,i),(i,i),..,(i,i) = 1 - 2(n - l)e (5.5) 

is, apart from the identity, the only relevant operator in the O.P.E. of e(x')e(x). The 
enlarged algebra of the fields and e is given by eq.(3.8), where Df e and Dt, are now 
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the related structure constant of the WD^ model. We consider thus the action (3.3), 
where in this case describes a single WD^ model: the correspondent /5-functions 
take the form (3.9) with the values of the structure constants (up to the first order in e) 
calculated in appendix 1: 



(5.6) 



^n(2n - 1) 

In the case N = 0, in addition to the unstable trivial fixed point (0, 0), there is one stable 
fixed point 

g, = 0, A, = 2y / n(2n- l)e (5.7) 

For N > 2, there are two other supplementary fixed points, both node points of the RG 
flow: 



A« = v /M2^)^l-(iV-2)^^- T 



n(n — 1) 



+ n(n- 1)N 2 



= +2(2n-l)e, 



n(n — 1) 



\ N - 1 + n(n - 1)N 2 



(5.8) 



A?) 



,(2) 



= y/n(2n-l)e ^1 + (iV - 2)^ ^-y 



n(n — 1) 



+ n(n-l)N 2 



= -2(2n-l)e„ 



n(n — 1) 



(5.9) 



'\| N - 1 + n(n - 1)N 2 

In order to compute the central charges at the fixed points for the disorder problem and 
for the coupled models, we can directly apply the result (4.9), provided the new values of 
the structure constants (5.6) and of the dimension of operators (5.4) and (5.5) are taken 
into account. For the quenched problem, the central charge c^] at the stable fixed point 
(5.7) is 



c(0,2 x /n(2n-l)e) 

,/;, - lim V N - = c(WD^)-An(n-l)(2n-l)e 3 ^c(WDt 1) ) (5.10) 



» 



The correspondent field theory is described by the WD^ ^model, as found already in 
[1]. In the coupled model, the central charge c^ upl at the same fixed point is c^ upL = 
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Nc(WD% ^) and the associated critical theory is described by N decoupled WD^ ^ 
models. 

At the new fixed points (|5.8|) and (5.9) the central charges cfy turn out to be: 



cf] = Nc(WD ( p } ) - 2Nn{n - l)(2n - 1) 



>(p) 



1 T (N - 2) 



n{n — 1) 



\ N - 1 + n(n - 1)N 2 



e 3 (5.11) 



As the RG-flow behavior is the same for all n > 3, the results we have shown in section 
4.1 for the disordered WD^ model and for the system of iV coupled WD^ models are 
valid for the whole series of WDffl conformal theories. 

Finally we give the direct generalization of eq.(4.6): 



A (i),(2) = Ae + (2 n - 1) J^l ± N, 
A (i),(2) = A £ + (2n-l) [ 1^N A 



n{n — 1) 



\ iV — 1 + n(n — 1)A^ 2 



n{n — 1) 



\ N - 1 + n{n - 1)N 2 



(5.12) 



6 Conclusions 



As has been stressed throughout the paper, for every small value of e: 

1 



p+ 1 



p > 1 



(6.1) 



we have a unitary model, WD^\ associated to it. When a set of iV > 2 of such models 
are coupled and brought, by fine tuning of the couplings g and A, to one of the two 
newly found fixed points, (5.8) or (5.9), the corresponding critical theory should also 
be described by a unitary conformal field theory. As we have a fixed point for every p, 
we should have a unitary series of new conformal theories, with p being a parameter, 
accumulating towards p = oo. Moreover, we should have a unitary series for every n, of 
WDM,n = 3,4,5,.... 

These new theories, presently unknown, have to incorporate into them the permuta- 
tional symmetry Sn, N = 2, 3, .... This symmetry has to be incorporated into the chiral 
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algebra of these theories. As the symmetry is discrete, the natural suggestion would be 
that it should be represented by parafermionic currents. 

For N > 3, the group Sn is non-Abelian. Looking at the expression of the central 
charge, eq.(5.11), we observe that the correction term which we have calculated is non- 
rational for N > 3. 

If we accept the idea of non-Abelian parafermionic conformal theories, we have to 
accept also that these theories possess an infinite series of unitary models, labeled by 
p, according to the arguments given above. And then the formula (5.11) indicates that 
the central charge for these unitary models takes non-rational values. This feature is 
unusual. 

Being more precise, the argument for rational or non-rational values of the central 
charge c, and the dimensions of the operators, could be given as follows. 

If one assumes that c, which is a function of p, takes rational values for all integer 
values of p, this then requires that the function c{p) should have a simple rational form, 
like c(p) = Q(p)/P(p), where Q(p),P(p) are two polynomials of p with rational valued 
coefficients. If the function c(p) of such a form is developped in a series of 1/p, the 
coefficients will all be rational. 

Still, it should be admitted that we don't know in fact if the perturbative expansion 
in e ~ 1/p, which we define in this paper, represents a convergent series. The series might 
also be asymptotic, as it is often the case in perturbative expansions in field theory. This 
then allows for two possibilities. 

The first possibility will be that the series is convergent, representing an analytic 
function c(p), like this is two case for the perturbed minimal models in [5]. 

In this case, according to the argument given above, one could actually judge on 
rationality- non rationality with the coefficients of the expansion. 

Saying it again, for having rational values of c(p), for all integer p, a complicated 
function c{p) will not be allowed. 

Then the coefficients of the expansion will also have to be rational. 
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The second possibility will be that the series in e is not convergent, that it is only 
asymptotic. In this case our arguments do not apply and the exact values of the central 
charge might well be rational, in spite of the irrational coefficients of the expansion. 

In view of this second possibility, our conclusions on non-rationality should not be 
viewed as definite. 

Construction of the corresponding exact conformal field theories represents a theo- 
retical challenge. 
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A Appendix 1 



A.l Relevant Operators in the O.P.E. of e{x')e(x) 

In section 3 we have explained how to deduce the O.P.E. of e(x')e(x) from the integral 
representation of the energy four-point correlation function: within the set of vectors 
{2/3 £ + g a e a } where 51 = 0,1,2, q 2 = 0,1, and q 3 = 0,1, we search for the 

ones which decompose into positive integer numbers of basic vectors of the 

representation lattice, i.e.: 

3 3 1 - ri 

2(3 £ + a_ ]T Qae a = a. ]T -—^u a (A.l) 

a=l a=l 1 

with n' a > 1. In this case the vertex operator V 2 ^ +a ^ q g belongs to the physical 
sector of the Kac table and the correspondent primary operator (f>( n ' i)( n ' 2 ,i)(n' 3 , i) appears 
in the energy algebra. Using eq.(2.19) and eq.(2.20), the condition (A.l) is equivalent to: 

n[ = 3 - 2q l + q 3 + q 2 
n' 2 = 1 - 2q 2 + qi 

n' 3 = l-2q 3 + qi (A.2) 

Testing all the possible set of values 52,93), it's easy to see that the system (A.2) 
admits only three solutions; they correspond to the identity operator (for qi = 2, g 3 = 
1, 52 — I), to the operator <p (for q 1 = 1, 53 = 0, q 2 = 0) and to the operator </>(3,i)(i,i)(i,i) 
(for 51 = 92 = 93 = 0). According to eq.(2.31) A( 3j1 )( 1j1 )( 1j1 ) = 2: the only relevant 
operators in the energy algebra are the identity and (/>. 

Now, in order to verify that the algebra formed by these two fields doesn't contain 
any other relevant operator, we have to check the O.P.E. of <f>(x')(/)(x). Following the 
same procedure used before, it's easy to verify that, apart from identity, the field is 
the only relevant operator in the 0-algebra, as already discussed in [1]. We consider 
in this case the integral representation of the four-point correlation function G'(x) = 
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(0(O)0(x)0(l)0(oo)): 

G\x) oc /../(^W^^^W^-MVf^OVf (^)^-(ei)^"(6)^"(^i)^(^)> 

(A.3) 

where the variables A*2, £1, £2, ^1, ^2 are integrated over the 2D plane. The charge 
neutrality is satisfied according to —2(3^ + 2a„(e 1 + e*2 + 63) = 0. As discussed above all 
the operators (ft (n'^i)^, 1)^,1) present in the O.P.E. of <fi(x')(f)(x) are given by the condition 

3 3 1 - n' 

2/^ + «_ (la^a = «- E ~ (A.4) 
a=l a=l 1 

where 91, q 2 , q 3 can assume the values 0, 1, 2; eq. (A.4) is equivalent to: 

n[ = 1 - 2qi + q 3 + q 2 
n' 2 = 3-2q 2 + q 1 

n' 3 = 3 + gi -2q 3 (A.5) 

By eq. (A.5) we have determined, apart from the identity, four primary operators: 0(i,i)(3,i)(3,i) , 
0(2,i)(3,i)(i,i), 0(2,i)(i,i)(3,i) and 0; the only field whose dimension is smaller than the unity 
is 0. 

A.2 Relevant Operators in the O.P.E. of e(x')e(x) for the WD$ 
model with n > 4 

As already discussed in section 5, the energy operator e of the WD^ model is associ- 
ated to the primary operator 0(2,i)(i,i)..(i,i)- We show below that in the correspondent 
algebra there is, apart from the identity, only one relevant primary operator, namely 
0(i,i)(2,i)(i,i)-(i,i) = 0- In the following we define for simplicity /3(2,i)(i,i)-(i,i) = fie, 
2a - /3(2,i)(i,i)...(i,i) = Pe, /3(i,i)(2,i)(i,i)-(i,i) = and 2a - /9(i,i)(2,i)(i,i)-(i,i) = 
The integral representation of the correspondent four-point correlation function G(x) = 
(e(0)e(x)e(l)e(oo)> is: 

G(x) oc / • • • / (^(0)^(^)^(1)^(00) Vf--- Vf V 2 -..-V 2 -...V-... V-) (A.6) 

mitimes m2times runtimes 
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where the screenings V a with a = 1, • • • , n are integrated over the 2D plane. Using the 
quadratic form (5.2) it's easy to see that the charge neutrality condition 

n 

-2/3 £ + a_£m a e a = (A.7) 

a 

imposes m a — 2 for a = 1, • • • , n — 2 and m n _i = m n = 1. For the same arguments 
discussed in the previous section, we search for the vectors ^(n'^Vj^iy-in'^i)) with n' a > 1 
(see 5.1) such that: 

n 

2(3 £ + a_ Vae a = ( A - 8 ) 

a=l 

where q a = 0, 1, 2 for a = 1, • • • , n — 2, q n _x = 0, 1 and q n = 0, 1. Using eq.(5.2), the 
condition (A. 8) is equivalent to the following system of equations: 

n\ — 3 - 2qi + q 2 

n' a = 1 - 2q a + g a _i + q a+1 2 < a < n - 2 

n 'n-2 = 1 ~ 2 <?n-2 + + + Qn 

n n _ x = 1 - 2g„_i + g„_ 2 

< = 1 - 2q n + q n _ 2 (A.9) 

Eq.(A.9) gives us all the possible primary operators 4>( n ' 1 ,i){n' 2 ,i)---{n' n ,i) present in the en- 
ergy algebra. Once again, testing all the possible set of values (?i, • • • , q n ) it's easy to 
verify that, apart from the identity, there are only two solutions: they correspond to the 
operator 0(3,i)(i,i)...(i,i) (for q a — 0), with dimension greater than unity, and to the oper- 
ator 0(i,i)(2,i)(i,i)-(i,i) = (for q± — 1 and q a = 0, a > 2), whose dimension is given by 
eq(5.5). In the operator product decomposition <f>(x')<f>(x), we have checked that the only 
relevant field is itself: by studying the correspondent four-point correlation function 
we obtain the following system: 

n[ — 1 - 2q 1 + q 2 

n' 2 = 3-2q 2 + q 1 + q 3 

n' i = l-2q i + qi_ x + q i+1 3 < i < n - 3 
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n' n _ 2 = 1 - 2g n _ 2 + q n _ 3 + q n _ x + q n 
rt n -\ = 1 - 2 <?n-l + Qn-2 

< = l-2g n + g n _ 2 (A. 10) 

with the integers qi, g n _i, g n = 0, 1, 2 and g a = 0, 1, 2, 3, 4 for a = 2, • • • , n — 2. According 
to eq.(A.lO), the field is the only operator with dimension smaller than unity which 
appears in the 0-algebra. 



B Computation of the structure constants Df e and 
for the WD$ models 

The structure constants Df e and Dt* are determined by two three-point correlation 
functions of the unperturbed WD^ model: 

Di = ( £ (O)5(l)0(oo)) 

^ = (0(O)0(l)0(oo)) (B.l) 

where corresponds to the primary operator 0(i,i)(2,i)(2,i) for the WD^ model and to the 
primary operator 0(i,i)(2,i)(i,i)--(i,i) for the WD^ model with n > 4. In order to compute 
their values using the Coulomb-Gas representation, we have to take into account that 
the vertex operators , V^_, and can acquire non-trivial normalization factors 
N £ , N g , and @] such that: 

m = N^(x) = N^V _(x) (B.2) 

After imposing the fields normalizations (5(1)5(0)) = N e N g = 1 and (0(1)0(0)) = 
N^N^ = 1, we need to compute four integrals: 

7 i = nMi)7y ■■J(V0,(O)V^(l)V^oo)^^-'-'^ : ^=N^ (B.3) 
* m * ' (1) *■ (1) +■ 

m\ times m n ' times 
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(2) 

n I \m 



m i 



n ( \m' 3 ' 



(2) t - (2) ,. 

m x times m n times 



i H I \IIL i> i> 



(4) 

n f \m 



m^ 3 ' times m^ times 



^^j--JW)VfiSW (B.6) 



m^ 1 times mj, 4 ' times 



where the integration over the 2D plane of the screenings, whose number is fixed by the 
charge neutrality condition, is intended. In eq. (B.6) the vertex corresponds to the 
identity operator / and the relative normalization constant Na ^ — Nj — 1. According 
to the general form of the correlation function of iV vertex operator 

w&) ■ • • *w&)> = n is - &r~ 4 > (b- 7 ) 

and using the formula 

C(a,f3) = /dWlC-l^ 

T(l + a)r(l + (3)T(-1 -a- /3) 



it- 



T(-a)T(-(3)T(2 + a + (3) 
and 



(B.8) 



K(a, P ) = /rf 2 ^ 2 ciene-i| 2a icric-i| 2a ic-^i 



ip 



2 r(2 P )r(i - P ) JL r 2 (i + a + ^)r(-i - 2a - (i + Op) , . 
71 r( P )r(i - 2 P ) H p(-a - ^p)r(2 + 2a + (i + i) P ) 1 • > 

where T(x) is the Gamma-function, the integrals (B.3)-(B.6) have been determined. 
Eq.(B.8) and eq.(B.9) represent special cases of a family of integrals calculated in ||. 

In terms of these integrals, the structure constants are : 



DZ, = I: 



I 



2 
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hh 



= (B.10) 
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B.l Computation of I\ 

The charge neutrality condition for I x is: 

-/^ + a_X>«e a = (B.ll) 

a=l 

It is satisfied for raj 1 ^ = = = 1, =2, a = 2, • • • , n — 2. Using eq.(B.7) 

and taking into account all the scalar products, easily computed from the quadratic form 
(5.2), the integral I\ takes the form: 

n-2 2 



/i = -2 3 ^/---/^inn^ ) ^n-Ax 



a=2 k=l 

2 n-3 2 

xi£ir 2a - n id fc) -eir 2Q -id fe) - i\- 2a - n n ief -&r 2 ^ x 



k=l a=2 k,l=l 

n-2 2 



x n itf > - e\ 4al n itfi - ^-ii- 2 - 2 - \& - u- 2a - ^ 

a=2 k=l 



where a 2 _ — p/(p + 1) ~ 1 — e. Integrating over the variables £ n , and £ n _i we obtain by 
eq. (B.8): 

n-2 2 



a=l k=l 

n-3 2 



xi6r 2a - n i# - 6r 2a -id fe) - ir 2Q - n n - diir 2a - x 

k=l a=2 k,l=l 



n— 3 

X 

a=2 



t — o 

U\e ) -er a -\&-&\ 4 - 4a - (B.i3) 



Then using eq. (B.9) we integrate in order over the couple of variables ^-21 " ' ' > £2^ an d 
over the variable £1; the result is: 

n— 3 
i=l 

B.2 Computation of I2 

According to : 

-2^ e + ^ + a_eJ = (B.15) 
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we have -* = 1 and — for a — 2, • • • , n; the integral (B.4) reads: 

h = -J d%(V^0)V^(l)V fr _(oc)V 1 -^ 1 )) 

= -/^iier 2Q -|6-ir 2a2 - 

= -C(-a 2 _,-a 2 _) (B.16) 
B.3 Computation of 7 3 

The integral (B.5) must satisfy the charge neutrality condition (B.ll); it takes the form: 

n-2 2 



-2 3 -« / . . . / d % n n d 2 e ] d%-id% 

J J a=2 k=l 

n - 6r 2a -id fe) - ir 2a2 -id fe) i- 2a - n n ie } - eSii- 2Q - x 



2 n-3 2 

X 

fc=l a=2 fcJ=l 



x ff \e ] - e ] r- n \& - m-**- \& - u- 2a - m?) 

a=2 k=l 

Following the same order of integrations used in the computation of the integral (B.3), 
with the difference that we integrate first over the variable £i and then over the couple 
of variables we obtain: 

I 3 = _2 3 -"C 3 (-a 2 _, -a 2 _)K(-a 2 _, —-^ - (n - 2)a 2 _) JJ K{-a 2 _,i{\ - a 2 _)) (B.18) 

2 i=\ 

B.4 Computation of I4 

According to the charge neutrality condition: 

n 

-2& + X>l 4 >e a = (B.19) 

a=l 

we have — 2 for a — 1, • • • , n — 2, m^li = = 1. The integral (B.6) has the form: 

h = 2 2 - n [■■■ fnn^'^-i^x 

J J a=l k=l 



2 n-3 2 

X 

fe=l a=l fc i=l 



n - i]- 2 «-id fc) r 2 « 2 - n n k?> -eSii- 2 -- x 
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n-2 2 

X 

a=l k=l 



n \& ] - & 4a - n iefi - a-ii- 2a - 1^2 - u~ 2a - (B.20) 



Integrating over the variables £ n , £ n _i, £„,_2? " * ' ? £1 j the integral (B.6) assumes the value: 

J 4 = C 2 (-a 2 , -a!) J] K(-a?_,i(l - a 2 )) (B.21) 

i=i 

Developing in e the integrals we have computed and using eq. (B.10), the eq.(5.6) are 
obtained. 



C RG Equations and Renormalization of energy op- 
erators 

Using the O.P.E.(3.8) and the dimensions of the energy and fields (eq.(5.4) and (5.5)), 
the operator algebra of the perturbing terms is: 

N N N 

E (e a e b ) (x) E (e c e d ) (y) - A(N - 2)\x - y\'^'^ £ (e a e b ) (y) + 

N 

+4(iV - i)i£|x - y\- 2+4n E ^a(y) + • • • 

a=l 

N N N 

E to) (*) E ^a(y) - 2^ E |x - yr^C-D E to) (y) + • • • 

a=l a^6 



N N N 

$>.(*) X> B (v) - D%\x- y \- 2 +^Y,My) + --- (c.i; 

a=l b=l a=l 

where we have omitted the irrelevant terms. By eq.(C.l) the 1-loop RG-equations can 



be easily obtained (see for example or [T(| ) 



& = 2(2n-l)eg-4n(N-2)g 2 -47tDt £ \g 

h = 2(n-l)eX-7rD%X 2 -27r(N-l)Dig 2 (C.2) 

With the redefinitions g — > g"/(47r) and A — > \/(2n), we find the eq.(3.9). 
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Similarly the renormalized energy operators e' c {x) (3.11) can be computed at the first 
order || from the following operator product decomposition 

N N 

£ (e a e b ) (x)e c (y) -> 2\x -y\-^^ 1 ^e a (y) + • 

a^b a^c 
N 

E Mx)e c (y) - I>£ \x - yr 2+2(n - 1)e e c (y) + • • • (C.3) 
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Figure 3: RG-trajectories of N = 3 coupled WDf' models. 
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